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Six years of XMM-Newton observations of NGC 1313 X-l and X-2 
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We present a systematic analysis of the X-ray spectra of NGC 1313 X-l and NGC 1313 X-2, using six years of XMM- 
Newton observations (17 observations). We fitted the continuum with a Comptonization model plus a multicolor blackbody 
disc, that describes the effects of an accretion disc plus an optically thick corona. We checked the consistency of this 
spectral model on the basis of the variability patterns of its spectral parameters. We found that the two sources show 
different spectral states. We tentatively interpret the observed behaviour of NGC 1313 X-l within the framework of super- 
Eddington accretion and that of NGC 1313 X-2 within the framework of near Eddington accretion. We also attempted to 
determine the chemical abundances in the local environment of NGC 1313 X-l and X-2, analyzing the RGS spectra of the 
longest observation available (122 ksec). The results appear to indicate solar metallicity for the local environment of NGC 
1 3 1 3 X- 1 and sub-solar metallicity for NGC 1313 X-2. 
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1 Introduction 



The barred spiral galaxy NGC 1313 hosts three Ultra lumi- 
nous X-ray sources (ULXs), one of which (NGC 1313 X-3) 
is a known supernova (SN 1978K) interacting with the cir- 
cumstellar medium. The other two ULXs, NGC 1313 X-l 
and X-2 (X-l and X-2 hereafter), are located in different 
positions, close to the nucleus X-l and at the outskirts of 
the host galaxy X-2. NGC 1313 has been observed several 
times by XMM-Newton and a sufficient number of X-ray 
spectra are now available to attempt a characterization of 
their spectral variability. 

The first analysis of the XMM-Newton spectra of X-l 
and X-2 was published by Miller et al. (2003). They fitted 
the spectra with an absorbed multicolor blackbody (MCD) 
plus a powerlaw model, finding that the temperature of the 
MCD component is much lower than that observed in Galac- 
tic black hole (BH) binaries. The high luminosity and cool 
disc component were interpreted as evidence for the ex- 
istence of intermediate-mass black hole (IMBH) of 100- 
1000M o (Miller et al. 2003, 2004). Feng & Kaaret (2005) 
also found that the early XMM-Newton spectra of X-l and 
X-2 were best fitted with a power-law plus multicolor disc 
blackbody model, commonly used to describe the spectra 
of accreting BHs. Later, Feng & Kaaret (2006) fitted a se- 
quence of 12 XMM-Newton observations with the same spec- 
tral model and found an anti-correlation between the disc 
luminosity and inner disc temperature. They concluded that 
the soft component does not originate in a standard accre- 
tion disc. Also Goncalves & Soria(2006) questioned the ro- 
bustness of the cool disc spectral component and, therefore, 
its use for estimating the BH mass. The optical and X-ray 



variability of X-2 was also investigated by Mucciarelli et al. 
(2007) to costrain the properties of the system. 

Finally, in an independent analysis, Stobbart et al. (2006) 
found that a possible curvature above 2-3 keV is present in 
high counting statistics spectra. More recently, adopting a 
more physically consistent Comptonization model and us- 
ing very long exposures, Gladstone et al. (2009) showed 
that many ULXs, including X-l and X-2, display distinct 
spectral curvature above 2 keV. This is interpreted as caused 
by an optically thick corona that hides the inner part of the 
accretion disc, in what is likely an extreme form of the so- 
called very high state of Galactic BH candidates (Done & 
Kubota 2006). 

In this work we try to characterize the spectral variabil- 
ity of X-l and X-2. Recently, a similar analysis was per- 
formed by Vierdayanti et al. (2010) on XMM and Swift data 
ofHoIXX-1. 



2 Data reduction 

We re-analyzed all the available XMM-Newton spectra of 
the ULXs hosted in the spiral galaxy NGC 1313, X-l and 
X-2, with homogeneous criteria. The 17 observations span 
a time interval of six years, from 17 October 2000 to 16 Oc- 
tober 2006. Three observations were excluded because of 
high flare contaminations. Data were reduced using SAS v. 
9.0.0. EPIC-MOS and EPIC-pn spectra were extracted se- 
lecting the good time intervals without background contam- 
ination and with a background count rate not higher than 
0.45 count s _1 . We used a 35" and 30" circular extraction 
region for X-l and X-2, respectively, and a 65" circular re- 
gion for the background. 
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RGS spectra were extracted using the rgsproc task with 
the option spectrumbinning=lambda. In this way it turns out 
to be possible to combine the spectra of different observa- 
tions. The EPIC spectra were grouped with a minimum of 
25 counts per channel, while RGS spectra with 50 counts 
per channel. 

All the spectral fits were performed using XSPEC v. 
12.5. 1 . To improve the counting statistics, we fitted the EPIC- 
pn and EPIC-MOS spectra together. EPIC spectral fits were 
performed in the 0.3-10.0 keV energy range, while RGS 
fits were limited to the 0.4-2.0 keV energy range. Follow- 
ing Gladstone et al. (2009), we modelled thermal Comp- 
tonization adoting the comptt model (Titarchuck 1994) or 
the eqpair model (Coppi 1999). When necessary, we added 
a multicolor black body disc model to the spectral fit (diskbb 
in XSPEC). The comptt model is an analytic model de- 
scribing Comptonization of soft photons, with a Wien in- 
put spectrum, in a hot plasma. The eqpair model allows 
for a 'hybrid' (thermal and non-thermal) plasma and com- 
putes the resulting comptonized spectrum without assum- 
ing that the electrons are non relativistic. The seed pho- 
tons may have a disc or blackbody spectral distribution. For 
ULXs, non-thermal processes are not likely to be impor- 
tant, so we decided to adopt a simplified version of eqpair, 
dubbed eqtherm, that neglects them. 

3 Results 

We analyzed 14 out of the 17 XMM-Newton observations 
of X-l and X-2 adopting an absorbed Comptonization plus 
multicolor blackbody disc components. The interstellar ab- 
sorption was modelled with the tbabs model in XSPEC. We 
fixed the Galaxy column density, along the line of sight, at 
0.39 • 10 21 cm" 2 (Dickey & Lockman 1990) and added a 
free absorption component to model local absorption near 
the source. We found that the X-l column density is sig- 
nificantly variable during time and goes from a minimum of 
1.1 -10 21 cm -2 to a maximum value of about 3.0- 10 21 cm -2 
for both comptt and eqtherm model. On the other hand, 12 
out of 14 observations of X-2 have values of the column 
density clustering around 1.4- 10 21 cm -2 . The other two ob- 
servations have N H = 0.3 • 10 21 cm -2 and N H = 0.5 • 10 21 
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Fig. 1 Optical depth versus temperature of the corona for 
the X-ray spectral fits of X-l (top) and X-2 (bottom) with a 
comptt+diskbb model. 



lies in the range 0.25-0.38 keV. When the soft component is 
present, the corona temperature is low (^ 2 kev) and the op- 
tical depth is high (r <: 7). When the soft component is not 
needed, 3 keV and r ~ 5. The spectra in these two 

different states are shown in Figured The only 'anomalous' 
observation without a soft component located in the region 
of high r-low kT cor in Figure \T\ (bottom) is one of the two 
with low Nh mentioned above. However, fixing the column 



cm l . If we fix the column density of these two observa- density equal t0 l A . 1Q 2i cm -2 and adding a soft compo _ 



tions equal to 1.4 • 10 21 cm 2 , the spectral parameters re 
main essentially unchanged. 

The best fitting parameters of the corona for X-l and X- 
2, obtained modelling Comptonization with comptt, are re- 
ported in FigureQ] All errors are at the 90% confidence level 
for one interesting parameter. The position of the observa- 
tions of X-l on the kT cor -r plot does not appear to show 
a clear correlation with the presence or absence of a soft 
component. The temperature of the soft component is in the 
range 0.21-0.35 keV. On the other hand X-2 appears to show 
a well defined behaviour in the kT cor -r plane that, apart 
from a single observation, is correlated with the presence 
of a diskbb component. The temperature of this component 



nent, the fit of this observation is statistically equivalent to 
that obtained applying only the comptt model (with fixed 
Nh), and the normalization of the soft component is com- 
parable to that obtained from the other spectra that show the 
presence of a disc. Therefore, it is unclear if this observa- 
tion is really anomalous or if there is not sufficient statistics 
a low energies to pinpoint the soft component. 

Letting aside this observation, when the disc component 
is present the total luminosity of X-2 may vary significantly, 
reaching ~ 8 • 10 39 erg s _1 . However, at low luminosity 
the source appears to show spectra with and without a soft 
component, likely reflecting the lack of one-to-one COrre- 
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spondence between spectral shape and luminosity observed 
in Galactic BH binaries. 




Fig.2 X-ray spectra of X-2 for a low luminosity state 
without disc (red) and a high luminosity state with disc 
(black). 

Figure [3] shows the light curves of X-l and X-2 com- 
puted from the best fitting diskbb+comptt models (a dis- 
tance of 3.7 Mpc was assumed; Tully 1988). On average, 
X-l has a luminosity (~ 10 40 erg s _1 ), higher than that of 
X-2 (~ 5 • 10 39 erg s _1 ). Variability of a factor ~ 3 and ~ 5 
is observed for X-l and X-2, respectively. 

We further investigated the behaviour of the soft com- 
ponent in both sources. X-l shows the disc component only 
in 5 observations while X-2 in 9 observations. In X-l the 
luminosity of the soft component is ~ 10 39 erg s _1 and it 
represents a significant part of the total flux. There is essen- 
tially no correlation between the luminosity of the disc com- 
ponent and the inner disc temperature. The same holds true 
also for the total luminosity. On the other hand, for X-2 the 
luminosity of the soft component is only a small fraction of 
the total luminosity and the disc luminosity (Ldi SC ) appears 
to show a power-law correlation with the inner temperature 
(Tin)'- L disc ex T 5 3±11 (see Figure In this fit we ne- 
glected the observation with kT cor = 1.7 keV and r = 10 
because it has a very low statistics and the disc parameters 
inferred from it are not well constrained. 

The results obtained modelling Comptonization with the 
eqtherm model are similar to those obtained with the comptt 
model for both sources. The soft component is needed in 5 
spectra for X-l and 7 for X-2. Also in this case X-2 dis- 
plays a dependence of the disc luminosity on the disc inner 
temperature, that follows the relation Ldisc T 21 ± - 5 . 

We analyzed also the RGS data of X-l and X-2 in an at- 
tempt to use them for determining the chemical abundances 
in the local source environment. As the sources are faint, 
only the last observation (the longest with 122 ks) reached 
a reasonable statistic for the brightest ULX (X-l), although 
the RGS net count rate is quite low (1.5 • 10 -2 count s _1 ). 
The continuum was fitted with the values inferred from the 
EPIC fit. We identify in the spectrum 2 lines in absorption, 



associated to the O I (535 keV) and Fe I (0.709 keV). There 
may be also two other lines in emission at the characteristic 
energy of O VIII K a (0.653 KeV) and Si K a (1.748 KeV), 
but their significance is very low. 

The lines do not allow to trace the chemical compo- 
sition. Therefore we substituted the tbabs model with tb- 
varabs, that allows to vary the chemical abundances and 
grain composition. We determined the abundances fitting 
the EPIC continuum and then fixing the abundance of a cer- 
tain element to zero. The RGS spectrum was then fitted with 
this model plus an absorption edge. The parameters of the 
edge are then used to compute the abundance. We did this 
with the Oxygen line at 0.538 keV and found an abundance 
consistent with solar. This result is confirmed repeating the 
analysis with the EPIC spectra, using the same method. 

We performed also a similar abundance analysis on the 
EPIC spectra of the last observation of X-2, finding a sub- 
solar abundance (~ 0.5 ± 0.1 Z ) for Oxygen while the 
Iron abundance is consistent with zero. This is in agreement 
with the abundances inferred from the HII regions in NGC 
1313 (e.g. Pilyugin 2001; Hadfield & Crowther 2007) and 
from the bubble nebula around X-2 (Ripamonti et al., these 
Proceedings). 

4 Discussion 

The analysis of all the XMM-Newton observations shows 
that the spectra of X-l and X-2 can be well reproduced by 
a Comptonization model plus a soft component. X-2 ap- 
pears to show a well defined spectral behaviour in the op- 
tical depth (r) versus corona temperature (T cor ) plane that, 
apart possibly for one observation, correlates with the pres- 
ence of a soft component. 

As already suggested by Gladstone et al. (2009), it is 
likely that X-l and X-2 are in different regimes. X-l shows 
higher average isotropic luminosity (~ 10 40 erg s _1 ) and 
smaller variability during the six years of observation. From 
time to time it shows the presence of a soft component that 
emits a significant fraction of the total flux and does not cor- 
relate with the disc temperature. As Gladstone et al. (2009) 
proposed, X-l could be in the ultraluminous regime, accret- 
ing at super-Eddington rates and launching powerful winds 
from the disc/corona. The accretion disc may be there, but 
covered by the wind/corona, while the soft component may 
actually represent emission from the wind itself. 

On the other hand, X-2 could be in a different regime, 
in which the average accretion rate is at around the Edding- 
ton limit. The disk may be still partly visible. This hypoth- 
esis is enforced by the correlation that we found between 
the disc luminosity and inner temperature, consistent with 
that of a standard disc (L oc T 4 ) when adopting the comptt 
model. The correlation disappears using a simple power-law 
to describe the spectrum of the corona, as found by Feng 
& Kaaret (2006), because it does not take into account the 
spectral curvature at high energies. The corona may become 
more expanded in the low-luminosity no-disc state, when it 
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Top: X-ray light curve of X-l (left) and X-2 (right). Unabsorbed luminosities evaluated in the 0.3-10 keV energy 



covers the entire disc and no disc component is visible in the 
spectrum. Clearly, as the corona is always optically thick, it 
is not possible to use the disc parameters for estimating the 
BH mass. We found also that for X-2 the corona temper- 
ature decreases and the optical depth increases as the total 
luminosity goes up, in agreement with the recent analysis of 
Vierdayanti et al. (2010) for Ho IX X-l. So X-2 could be in 
a state similar to that of Ho IX X-l, with the corona becom- 
ing progressively mass loaded as the luminosity increases 
(Gladstone et al. 2009). 

A couple of ULXs show periodic intensity variations in 
X-rays which are considered as signatures of the orbital pe- 
riod (Kaaret et al. 2006ab; Kaaret & Feng 2007; Strohmayer 
2009). X-ray spectra seem to change regularly with orbital 
phase. It is possible that these phase related variations may 
affect any result obtained from snapshot observations such 
as presented here. This can be addressed only through a ded- 
icated X-ray monitoring programme for X-l and X-2. 

We used the RGS high spectral resolution to attempt an 
estimate of the metallicities of the local environments of X- 
1 and X-2. Because of the low signal-to-noise ratio, only the 
last, longest observation of X-l could be used. The analysis 
was performed also on the EPIC spectrum. The metallicity 
in the X-l environment is consistent with being solar and 
appears to be higher than the average metallicity of NGC 
1313. The last EPIC spectrum of X-2 was also analyzed 
with the same method and suggests sub- solar metallicity. At 
present we are working at stacking together the RGS spectra 
of the longest XMM-Newton observations (performed in the 
same state for X-2), trying to improve the signal-to-noise. 
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